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We model diffraction errors found when us- 
ing toothed apertures [L. P. Boivin, Reduc- 
tion of diffraction errors in radiometry by 
means of toothed apertures, Appl. Opt. 17, 
3323-3328 (1978)]. Using toothed (cf. 
circular) apertures minimizes diffraction by 
inducing destructive interference within the 
diffracted signal. Since diffraction effects 
can be quite complicated, their over-all 
reduction may help limit uncertainties in, 
say calibrations. Our analysis yields three 
principles to guide design of nonlimiting 
(baffle) apertures which minimize diffrac- 



tion. We performed detailed diffraction 
calculations within scalar (Kirchoff) dif- 
fraction theory, using parallel-computing 
resources at the National Institute of 
Standards and Technology. 
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1. Introduction 

Limiting or nonlimiting apertures (the latter of which 
are also called baffles) and occulting disks are used in 
applications ranging from radiometric calibration [1^] 
to solar coronagraphy [5]. However, the utility of these 
devices is hampered by imperfect knowledge of diffrac- 
tion of radiation at their edges. In radiometry, for in- 
stance, diffraction leads to deviation from geometrical 
optics in the total radiation, from a given source, inci- 
dent on a detector. Here, we consider the technique of 
toothing aperture edges to reduce diffraction effects, 
being motivated by Boivin's demonstration of the effi- 
cacy of this approach [4]. We model Boivin's experi- 
mental results, and we formulate and test three princi- 
ples to guide design of apertures which lead to minimal 
diffraction effects. Specifically, we discuss only the 
effects of teeth on radiation detected . 

We limit this work to optics having broad-band (e.g., 
thermal) sources and fully illuminated detectors sepa- 
rated by screens with nonlimiting, (toothed) circular 



apertures. Diffraction of radiation from broad-band 
sources is manifested by the total detected flux differing 
from that predicted by geometrical optics. One desires 
to know the ratio of actual flux to "geometrical" flux. 
This ratio is often called "F2" for optics like those stud- 
ied here, and "F/' for certain, analogous optics having 
limiting apertures [1^]. As one is concerned with the 
difference between these ratios and unity, a ratio (Fi or 
F2) may be reexpressed as 1 -1- {s). This {s) is found by 
appropriate integration, over wavelength, A, of the rela- 
tive difference in flux incident for each A, ^(A), which 
could be positive or negative. 

In the "F2" case, Boivin found large effects on {s) 
from both depth and frequency of teeth, and noted how 
Huygens' principle suggested only a tooth's aspect ratio 
would affect diffraction. Better analysis of diffraction 
suggested a minimum necessary tooth depth for 
a significant reduction in {s). In general, sufficiently 
deep teeth reduce edge-diffraction effects through path- 
length-related phase cancellations between diffracted 
rays created at different points on an aperture 
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perimeter. This destructive phase interference is the key 
mechanism responsible for Boivin's observation of a 
reduction in diffraction effects. However, detailed calcu- 
lations to model the measured, reduced {e)'s were not 
carried out in Boivin's presentation. 

This paper models the diffraction effects measured by 
Boivin. The boundary-diffraction-wave formulation [6] 
helps clarify the effects of depth and frequency of teeth 
on diffraction. This leads to guiding principles for de- 
sign of toothed apertures optimal for minimizing dif- 
fraction. We test these principles in our own diffraction 
calculations. Our objectives include understanding the 
diffraction by toothed apertures, and designing such 
apertures to control diffraction. This work also demon- 
strates the feasibility of computing diffraction effects 
for complex optics, e.g., irregular apertures. We use a 
parallel-processing implementation of computer pro- 
grams which model diffraction. Numerical uncertainties 
in our results are controlled, isolating numerical approx- 
imations to the scalar Kirchoff theory (which we use) 
from physical approximations of that theory (or, in the 
case of comparison to measured diffraction effects, ex- 
perimental errors). In the Kirchoff approach, one evalu- 
ates the value of a scalar radiation field behind a screen 
using the Green's Function and the value of the undif- 
fracted radiation field incident on an aperture. 

Below, we first discuss diffraction theory and compu- 
tational issues pertinent to this work. We next present 
diffraction calculations modeling Boivin's {e)'s for cir- 
cular and toothed apertures in 25 optics. Then we iden- 
tify and apply guiding principles for design of apertures 
exhibiting minimal diffraction, testing aperture designs 
within Kirchoff theory by further calculations. We close 
with some conclusions. 

2. Diffraction Theory, Computational 
Issues 

Consider (cf. Fig. la) an optic consisting of a circular, 
extended source (radius p) and fully illuminated, circu- 
lar detector (radius r) placed on opposite sides of a 
screen with an aperture (nominal radius R). Source- 
screen and screen-detector distances are respectively a 
and b. Source, aperture and detector areas lie within 
parallel planes, and the centers of the areas are colinear, 
defining an optical axis. R is a "nominal" radius only, 
because we might consider a toothed aperture. The 
perimeters of toothed apertures are defined as follows. 
One begins with a circular aperture with radius R , and 
cuts A^ teeth by forming a perimeter consisting of 2N 
straight line segments. At intervals separated by angle 
(j) = 3607(2A^), which is subtended by half of a tooth, 
the aperture radius alternates between the values, R and 
R-\- A, where A is the tooth depth (cf. Fig. lb). 



Aperture 
Source j Detector 

I^ llL. 



-PP 

< a ^ b > 



Optic Axis 



Fig. 1(a). Schematic diagram of F2 optic: extended source, screen 
with aperture, and detector; various physical dimensions are indicated. 




angle, <^ 



Fig. 1(b). Section of toothed- aperture perimeter. 




^^^^^^M. 



Fig. 1(c). Geometry involved in the boundary- diffraction- wave 
formulation. Note points on source (rs), on detector (ro), on aperture 
perimeter (rA), and related vectors, s, t, and t\ 



Screen 




Aperture 



\ Trace of simple teeth 

Trace of circular arc 



Fig. 1(d). Section of novel, toothed- aperture perimeter. 
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To evaluate (s), one could perform extensive integra- 
tion over wavelength and areas of the source, aperture 
and detector. Use of an "effective-wavelength" [1-3] 
eliminates wavelength integration, and the boundary- 
diffraction-wave formulation of Kirchoff theory re- 
places double integration over aperture area with single 
integration over its perimeter. Symmetry simplifies inte- 
gration over either the source or detector. Consider such 
integration in terms of polar coordinates. For instance, 
treatment of diffraction by a circular aperture requires 
no angular integration over the source, and requires an- 
gular integration over only half of the detector. Treat- 
ment of diffraction by a regularly toothed aperture re- 
quires angular integration over only a narrow angular 
wedge (pie slice) of the source, and integration over only 
half (or all) of the detector. Integration over the source 
can be avoided in cases where assuming an axial, point- 
source is a valid approximation. Then one also needs to 
integrate flux over only a small part of the detector. We 
were not able to establish a predictive criterion for deter- 
mining whether this simplification affected results. 

Below, we discuss the formula used to describe 
monochromatic flux incident on the detector, whence a 
formula for £(\). We also discuss the boundary-diffrac- 
tion-wave formulation, and the effective-wavelength 
approximation used in this work, touching finally on 
considerations regarding the fineness of integration over 
points on the source, aperture, and detector in the 
present optics. 

We use scalar diffraction theory, so the radiation field, 
^, obeys the Helmholtz equation in free space. 



Here, r^ samples all the aperture area. The derivatives 
are taken normal to the screen plane in the direction 
towards the detector; ^i is the incident radiation field, 
not yet specified. The parameter jS plays the following 
role: Having p = yields the Kirchoff theory; having 
jS be +1 or -1 respectively yields extensions of Kirchoff 
theory which more consistently satisfy '^= or d'^ldn 
= boundary conditions on the dark side of the screen; 
cf. Jackson [7] for further discussion of these points. 

For a scalar radiation field, intensity is related to a 
radiation current density. 



J(ru) = ^ [^*(rD) Vnru) - ^(ro) V^^ir^)] . (4) 



However, for our optics, this is approximately propor- 
tional to I ^(ro)!^. The relative error in the computed 
intensity, if the latter is estimated using such an approx- 
imate value, is \/kb, or typically 10"^, because kb is 
typically 10^. Therefore, we simply use the square of the 
radiation field when computing intensity. 

If one has an incoherent, extended, broad-band source 
such as a lamp, detected radiation at each wavelength is 
often expressed as the sum of fluxes from many 
monochromatic, mutually incoherent point sources, and 
we compute our total flux as such a sum. Consider a 
monochromatic, point source located at point r^ on the 
source. Its radiation field as determined from geometri- 
cal optics, which we denote as '^g, is given by the 
following rule. In the parts of space illuminated by such 
a source, one has 



[V + k^] ^(r) = 0. 



(1) 



'^G(rD)| ilium = 



^ exp[/^|rD-rs|] 



(5) 



Here k is the wave number, i.e., Itt/X . More specifically, 
we use Kirchoff theory, which solves for ^on the detec- 
tor's side of the screen using the Green's Function 
approach. The Green's Function is given by 



, exp [ ik\r-r'\] 
^ ' ■' 4 Tr\r-r'\ 



(2) 



However, in the geometrical shadow of a screen, one has 



V^G(rD)|shad = 0. 



(6) 



In the boundary-diffraction-wave approach, one 
rewrites the full ^ arising from this point source, in- 
cluding diffraction effects, in this simplified form: 



The radiation field at a point Td, on the detector 
(Fig. Ic), is given as 



nrj,) = ^G(r J + 



/ 

aperture 



dVA 



nr,) = 



(13+1) Giro, Ta) — ^i(rA) + 



(13 - 1) ^(rA) ^ Giro, Ta) 



(3) 



f . [[ s Xt \ c\p[h 



kjs+t)] 



Amt 



_ n( ^Xf \ cxpUkjs + t')] '] 
'^\st'+s-t') 4mt' J ■ ^ '' 
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This is exactly equivalent to the previous formula. Here, 
s and t are respectively vectors from points rs and Td to 
Ta; t' is the vector, to point Ta, from the point obtained 
by reflection of Td through the screen plane. These vec- 
tors are illustrated in Fig. Ic. For geometries where s 
and t are nearly anti-parallel, the value of ^ is irrelevant. 
For our optics, this is true; expediency dictates that we 
set )S to zero, simplifying diffraction calculations. For 
later purposes, we make the abbreviation. 



nr) = ^G(r) + ^bC/-). 



(8) 



Furthermore, we define ^b(''d, /"s; A) as the value of 
^B (/"d) arising from a point source at rs emitting radia- 
tion at wavelength A; there is an analogous function, 
^gCz-d, A). 

Next, for a spatially incoherent, monochromatic, 
extended source, we may write 



1+^(A) 
J dVs J dVD|^G(/'D,/'s;A)+^B(/'D,/'s;A)|' 

source detector 



/dV.[/ 



dVD|^G(rD,rs; A)|^ 



(9) 



The integrand in the numerator may be rewritten as 

|^G(rD,rs;A)+^B(rD,rs; A)p 

= |^G(rD,rs; A)|^ + |^B(rD,rs; A)|^ 

+ ^*G(rD,rs;A)^B(rD,rs; A) 



+ ^G(rD,rs;A)^*B(rD,r5;A). 



(10) 



Compared to the first (purely geometrical) term, the 
second term varies roughly as A for optics studied here. 
For optics in which diffraction of detected radiation 
implies a substantial angular deflection of that radiation 
by the aperture, the third (interference) cross-term 
oscillates rapidly with A and can be largely self- 
cancelling for broad-band sources. This was found to be 
the case in the present work. 

Such cancellation allows us to use an effective-wave- 
length approximation, and the particular approximation 
used here — which differs slightly from earlier effective- 
wavelength approximations — provides a A -averaged 
value of ^(A), (s), as follows: 



jdA^(A)5(A)D(A) 

(s) = —r ^o«A)). 



I 



(11) 



dA5(A)D(A) 



For cases when ^(A) varies as A, valid in this work, the 
effective wavelength is 



(A): 



dAA5(A)D(A) 



(12) 



dA5(A)D(A) 



S{k) and D{k) are the source spectral density and 
detector responsivity, respectively; £q{\) is implied by 
Eq. (10) after dropping the cross-term. For p«a+b, 
the bracketed integral in the denominator in Eq. (10) 
depends weakly on rs. Ignoring this dependence leads to 
a relative error of [p/(a+Z?)]^, i.e., typically 10'^. There- 
fore, we ignore this dependence. When using Eq. (10) to 
evaluate {s), or to determine a diffraction pattern for 
other purposes, one must perform many independent 
evaluations of ^^(rD, rs; A) and ^o (/"d, ''d; A) for differ- 
ent values of rs and rD. Use of parallel computing tech- 
nology can easily accelerate calculations, and we have 
exploited this fact. 

Geometrical and other considerations indicate the re- 
quired detail in sampling rs, ro, and rA. Ideally, one 
would conduct only the coarsest sampling necessary. 
Numerical convergence is improved at the cost of fur- 
ther computation. From calculations which included ei- 
ther coarser or finer samplings than were used to obtain 
the results presented, we estimate uncertainties, arising 
from controlled errors in numerical integration, as fol- 
lows. For each {s), call its expanded uncertainty 6(^). 
Then 8{s)/(s) is around 0.03 for the results in Table 1 
(simulation of Boivin's experimental numbers), and 
0.05 for the results in Table 2 (simulation of a novel 
design of toothed aperture, to be discussed in a later 
section), except for results for extended sources, where 
it is 0.10. These estimates do not include biases related 
to the effective-wavelength approximation. 

We sampled ro at 0.01 mm radial and 0.25° angular 
intervals, respectively. We improved sampling of rs in- 
crementally, stopping when we could use numerical in- 
terpolation to estimate diffracted flux for rs everywhere 
on the source. Consequently, the radial coordinate of rs 
was sampled at intervals of 0.05 mm or 0.1 mm. We 
never found substantial dependence of the diffracted 
flux on the angular coordinate of rs. Symmetry was 
extensively exploited to accelerate calculations. 

Two factors assisted convergence of results with re- 
spect to detail of integration. First, the (^)'s were sums 
of strictly positive numbers. Second, distances between 
various optical components were large compared to 
components' dimensions, which were transverse to the 
optical axes. So phases of emitted or diffracted radiation 
exhibited only gradual spatial variations over the aper- 
ture and detector areas. This facilitated sampling rA at 
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intervals up to 200 times A . Integration about the aper- 
ture perimeter in Eq. (7) was usually assisted by two-, 
four-, or six-point Gaussian quadrature over each line- 



segment sample of the perimeter. Use of progressively 
high-ordered Gaussian quadrature permitted diminish- 
ing returns of acceleration of the integration. 



Table 1. Diffraction effects for optical configurations discussed in text 



Config. 


a 


h 


R 


A 


<^ 


r 


(s), theory 


(s), expt. 




(cm) 


(cm) 


(mm) 


(mm) 


C) 


(mm) 


(%) 


(%) 


1 


50 


50 


3.5 


smooth 




1.25 


0.682 


0.74(5) 


2 


50 


50 


3.5 


0.01 


3 


1.25 


0.642 


0.59(4) 


3 


50 


50 


3.5 


0.02 


3 


1.25 


0.541 


0.52(4) 


4 


50 


50 


3.5 


0.05 


3 


1.25 


0.106 


0.18(4) 


5 


50 


50 


3.5 


0.09 


3 


1.25 


0.023 


0.04(4) 


6 


50 


50 


7.5 


smooth 




1.25 


0.310 


0.36(4) 


7 


50 


50 


7.5 


0.01 


3 


1.25 


0.288 


0.18(4) 


8 


50 


50 


7.5 


0.02 


3 


1.25 


0.234 


0.15(4) 


9 


50 


50 


7.5 


0.05 


3 


1.25 


0.0246 


0.03(4) 


10 


50 


50 


7.5 


0.09 


3 


1.25 


0.0046 


0.01(4) 


11 


85 


40 


3.5 


smooth 




1.6 


0.464 


0.32(4) 


12 


85 


40 


3.5 


0.1 


9 


1.6 


0.374 


0.32(4) 


13 


85 


40 


3.5 


0.1 


5 


1.6 


0.149 


0.08(4) 


14 


85 


40 


3.5 


0.1 


3 


1.6 


0.026 


0.04(4) 


15 


85 


40 


3.5 


0.1 


2 


1.6 


0.017 


0.02(4) 


16 


85 


40 


3.5 


0.1 


1.5 


1.6 


0.014 


0.04(4) 


17 


50 


85 


3.5 


smooth 




1.25 


1.10 


1.05(5) 


18 


65 


70 


3.5 


smooth 




1.25 


0.95 




19 


80 


55 


3.5 


smooth 




1.25 


0.77 




20 


95 


40 


3.5 


smooth 




1.25 


0.57 




17+18 






3.5 


smooth 




1.25 


2.05 


2.06(5) 


17+18+19 






3.5 


smooth 




1.25 


2.81 


2.95(5) 


17+18+19+20 






3.5 


smooth 




1.25 


3.38 


3.52(5) 


21 


50 


85 


3.5 


0.1 


3 


1.25 


0.025 


0.00(5) 


22 


65 


70 


3.5 


0.1 


3 


1.25 


0.033 




23 


80 


55 


3.5 


0.1 


3 


1.25 


0.035 




24 


95 


40 


3.5 


0.1 


3 


1.25 


0.025 




21+22 






3.5 


0.1 


3 


1.25 


0.058 


0.06(5) 


21+22+23 






3.5 


0.1 


3 


1.25 


0.093 


0.08(5) 


21+22+23+24 






3.5 


0.1 


3 


1.25 


0.118 


0.13(5) 


25 


50 


50 


3.5 


1.4498 


45 


1.25 


0.13 


0.2 



Table 2. Diffraction for optics discussed in the text. 



No. of teeth 


Is principle 3 
applied? 


Source type 


A 
(H-m) 


(%) 


120 


no 


point 


0.58 


0.000 43 


120 


yes 


point 


0.58 


0.000 40 


240 


no 


point 


0.58 


0.000 13 


240 


yes 


point 


0.58 


0.000 23 


480 


no 


point 


0.58 


0.000 13 


480 


yes 


point 


0.58 


0.000 05 


960 


no 


point 


0.58 


0.000 13 


960 


yes 


point 


0.58 


0.000 000 43 


960 


yes 


extended 


0.71 


0.000 009 


960 


yes 


extended 


0.58 


0.000 008 


960 


yes 


extended 


0.49 


0.000 017 
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3. Modeling Measured Diffraction Effects 

Data modeled were obtained by Boivin [4] using a 
3000 K, 1 mm diameter, tungsten-radiation source and 
RCA 6217 photomultiplier as the detector/ The effec- 
tive detector area was controlled by 2.5 mm or 3.2 mm 
diameter proximity apertures. Boivin presented results 
for at least 27 optics involving one or more intervening 
apertures. The (A) for this source-detector combination 
was 0.58 |ULm. When modeling the data, we shall also 
analyze the effective-wavelength approximation. Be- 
sides that approximation, deviation of our results from 
experiment can be attributed to Kirchoff theory and/or 
nonideal experimental circumstances. Experimental dif- 
ficulties could include misalignments of optical compo- 
nents, even by 0.1 mm (based on our results), or could 
involve irregularities in aperture perimeters, an effect 
stressed in Ref. [4]. 

We computed theoretical (£^)'s for 25 of Boivin's 
optics. Results for combinations of apertures were ob- 
tained by adding (£^)'s computed for the individual aper- 
tures separately. (Effects of multiple diffraction by sev- 
eral apertures were assumed to be negligible.) Present 
results and results by Boivin are tabulated in Table 1. 
Also shown are parameters specifying optics. The (£^)'s 
are given as percentages. We also give the experimental 
uncertainties in (s), which are indicated in Ref. [4]. 
Config. 25 in Table 1 involved a 7 mm X 7 mm, square 
aperture. A 's reported in Ref. [4] were up to 0.003 mm 
different from those used here. This is a negligible 
effect. In Figs. 2 to 5, we plot results for Config. 1 to 
Config. 24 in the format of Figs. 3, 6, and 8 of Ref. [4]. 
In Figs. 2 and 3, we indicate traditional theoretical val- 
ues of (s) for a round aperture [1^], 






(13) 



as well as extrapolated values valid (within Kirchoff 
theory) for small A , 



<^)-<^) 



i^om 



(14) 



This latter formula follows from the observation that 
first-order effects of A , those because of a change in 
average aperture radius, are equivalent to a similar 
change in the radius of a circular aperture. 



^ Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 
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Fig. 2. The values of (s) for 7 mm diameter apertures with indicated 
tooth depth, cf. Configs. 1 to 5 in Table 1. Theoretical points are 
plotted with squares, experimental points are shown as lozenges. The 
dashed line shows the theoretical behavior of (s) in the limit of A 
approaching zero. 
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Fig. 3. The values of (e) for 15 mm diameter apertures with indi- 
cated tooth depth, cf. Configs. 6 to 10 in Table 1. Theoretical points 
are shown as squares, and experimental points are shown as lozenges. 
The dashed line shows the theoretical behavior of (e) in the limit of 
A approaching zero. 
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Fig. 4. The values of (s) for 7 mm diameter apertures with various 
numbers of teeth, versus the diameter of a "diffraction- free region," 
suggested by a geometrical model in Ref. [4]; cf. Configs. 11 to 16 in 
Table 1. Theoretical points are shown as squares, and experimental 
points are shown as lozenges. 
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Fig. 5. The values of (s) for combinations of apertures from 
Configs. 17 to 20 measured in Ref. [4], and for combinations of 
apertures from Configs. 21 to 24, also measured in Ref. [4]. The 
Configs. are described further in Table 1 . Theoretical points are shown 
as squares, and experimental points are shown as lozenges. 



4. Discussion 

We now present results assessing the effective-wave- 
length approximation, motivate several guiding princi- 
ples for the design of minimally diffracting apertures, 
and apply these principles in (simulations of) several "F2 
case" optics. 

4.1 Effective-Wavelength Approximation 

The present effective-wavelength approximation 
assumes that SoW is a local average of £r(A) with re- 
spect to A , and that this local average varies as the first 
power of A. We have computed SoW and s(\) for a 
range of A for Config. 8 in Table I, but with an axial 
point source, and we have computed SoW for a range of 
A for the same configuration, for a I mm diameter, 
incoherent source. These results are presented in Figs. 6 
and 7. Evidently, the effective-wavelength approxima- 
tion affects results by a small amount of {£:) for a point 
source, and to a lesser degree for an extended source. 
Selection of Config. 8 for the analysis was spurred by 
the large discrepancy between its theoretical and exper- 
imental {£^)'s. (Based on these tests of the effective- 
wavelength approximation, we conclude that such a dis- 
crepancy is not attributable to the effective-wavelength 
approximation.) 
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4.2. Principles for Minimizing Diffraction 

Consider an F2 optic with a circular, toothed aper- 
ture, where diffracted radiating reaching the detector 
would have experienced only a small angular deflection 
at the aperture perimeter, yet a sufficiently large deflec- 
tion to permit the present effective-wavelength approxi- 
mation. Assume that Kirchoff theory is a valid descrip- 
tion of diffraction in one's optic. In calculating ^b('*d), 
the vector defining an aperture-perimeter segment, d^A, 
has azimuthal and a radial components (cf. Fig. lb): 



dTA = aiU + a2i 



(15) 
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(We use cylindrical, polar coordinates, m, 0, and z, and 
still consider planar screens with apertures. Use of w for 
radial coordinate prevents confusing it with source 
radius, p.) 

When the aperture radius varies azimuthally — e.g., 
there are teeth — there will be path-length-induced phase 
differences between diffracted rays, so their contribu- 
tions to ^b(''d) by the azimuthal components of the 
dTA's will have large cancellations. Also, these contribu- 
tions will converge to some limit, if the depth of teeth is 
fixed, but the frequency of teeth is increased. 

Regarding analogous contributions to ^b(''d) by the 
M-components of the drA's, there should also be cancel- 
lations because of similar effects of phase. Also, cancel- 
lations arise from the alternating, inward/outward direc- 
tion of drA- Concomitantly, contributions to ^b(''d) from 
the M-components of the drA's will typically exhibit high 
azimuthal variations, and regularity of ^b requires con- 
tributions with high azimuthal variations to vary as high 
powers of u when u is small. Suppose there is a very 
high frequency of teeth on one ' s aperture. Then, in the 
case of a sufficiently remote and/or small source close to 
the optical axis, a small detector which is also close to 
the optical axis will receive minimal flux arising from 
the u -components of the d/^A's. We confirmed this effect 
numerically in our calculations. 

With the above considerations, two principles of de- 
signing baffle apertures that minimize diffraction are as 
follows: 

(1) Having deep teeth should help reduce diffracted 
flux in the central region, because of path- 
length-induced phase-cancellation effects. 

(2) Having a high frequency of teeth prevents 
contributions to the diffracted flux because of 
radial components of aperture-perimeter seg- 
ments. This permits one to consider diffracted 
radiation as the coherent superposition of dif- 
fracted radiation from several hypothetical, 
concentric, circular apertures with different 
radii. In such a picture, each radius is weighted 
according to the fraction of the actual aperture 
perimeter which is at such a radius (distance 
from the optical axis). 

Both of these principles are consistent with the results in 
Table 1. 

For high frequencies of teeth, the diffracted radiation 
field in the central region is related to a Fourier trans- 
form of the radial distribution of the aperture perimeter. 
Based on Parseval's theorem, relating the integral of the 
square of a function to the integral of the square of its 



Fourier transform, one anticipates that diffraction 
effects would be smaller for more deeply cut teeth. The 
diffracted flux is related to the integral, over radius, of 
the square of the fraction of the aperture perimeter at 
each value of radius. However, the integral over radius of 
the first power of that fraction must be one. 

One might wish to constrain tooth depth. Very deep 
teeth may be difficult to realize and, obviously, under- 
mine the functioning of a baffle. Given such a con- 
straint, we propose a third principle for design of mini- 
mally diffracting apertures: 

(3) Whereas diffracted rays from any single tooth 
might not be sufficiently weak, diffracted flux 
on a small, central detector may also be reduced 
by radially displacing half of the teeth by the 
distance 7rab/[{k)R(a-\-b)], (k) being 27tI{\). 
For diffracted rays reaching the detector, rays 
originating from one tooth will be similar to, 
except differing by a 1 80° phase-shift from, rays 
originating from an adjacent tooth (cf. Fig. Id). 

The 180° phase-shift depends on the wavelength and 
above, designed path-length difference. So this third 
principle works best when one's source and detector are 
smaller than an aperture, and the effective bandwidth of 
radiation transferred for the source-detector combina- 
tion is small compared to the radiation's effective central 
frequency. The 1 80° phase-shift being exploited would 
vary as the wave number of radiation from a broad-band 
source, somewhat degrading the level of destructive in- 
terference. The wavelength range having substantial 
interference would presumably encompass many oscil- 
lations with respect to A of the relative signs of the 
geometrical and boundary-diffraction waves. 

To illustrate these three principles further, we 
consider a set of optics like the others used in this work. 
We use a = b = 5^ mm, r= 1.25 mm, R = 7.5 mm, and 
either an extended source with p = 0.5 mm or an axial, 
point source. We vary the number of teeth, A^, from 120 
to 240, 480 and 960, and use several effective wave- 
lengths. Without application of the third principle, we 
set A = 0.2 mm (cf Fig. lb). Application of the third 
principle to aperture design is as follows (cf. Fig. Id). 
For the case of A^ teeth, divide an aperture into N/4 
equal angular wedges. Further, divide each wedge 
equally into 8 smaller wedges. The aperture perimeter 
contains one straight segment within each of these 
smallest wedges. At the nine, equally spaced angles, 
which collectively define both edges of all eight smallest 
wedges, the aperture radius has the following values: 
7.5 mm, 7.704 84 mm, 7.509 67 mm, 7.714 50 mm, 
7.519 34 mm, 7.714 50 mm, 7.509 67 mm. 
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7.704 84 mm, 7.5 mm. All apertures reflecting the third 
principle are designed for A = 0.58 |ULm. So we radially 
displace alternate pairs of teeth. In Table 2, we present 
results for (e)'s with different combinations of N, (A), 
source, and whether or not the third principle is applied. 
The results in both Tables demonstrate all three princi- 
ples, and the third principle is robust with respect to 
small changes in A. 

4.3 Miscellaneous Issues 

One should critically assess Kirchoff theory and the 
feasibility for manufacturing and using toothed aper- 
tures. This is particularly true for the apertures with the 
third principle applied, whose design relies on very 
small structure. One should also consider effects of edge 
roughness in real apertures. Within Kirchoff theory, suf- 
ficiently small features should not have large effects. 
However, since careful radial displacements of teeth by 
about 1 fxm dramatically influence diffraction, similar 
irregularities in the apertures of Ref. [4] might con- 
tribute to the discrepancy between our (e)'s and those 
reported in this work. 

5. Conclusions 

We have studied the approach of toothing apertures to 
reduce diffraction. We modeled diffraction effects found 
by Boivin [4] for 25 optics, each involving an extended, 
incoherent, broad-band source, one or more apertures, 
and a fully illuminated detector. Modeled diffraction 
effects agreed reasonably with those measured, and 
remaining discrepancies could arise from use of the 
Kirchoff approach or experimental nonidealities (aper- 
ture irregularities, in particular). We also discussed 
three principles which can guide design of apertures 
that minimize diffracted flux incident on a central 
detector. Successful application of these principles was 
demonstrated, not only in modeling diffraction effects 
measured by Boivin, but also in simulations of specially 
designed apertures. 

We carried out detailed diffraction calculations using 
the Kirchoff formulation of scalar diffraction theory, 
with the work being facilitated by the boundary-diffrac- 
tion-wave technique. We estimated the level of numeri- 
cal convergence of these calculations. An effective- 
wavelength approximation was used, and its validity was 
assessed. Calculations were accelerated by use of paral- 
lel processing, and acceleration of calculations should 
continue to scale well up to networks of large numbers 
of processors. 

Future research considered includes application of 
toothed apertures in radiometry, the role of aperture 
irregularities, further assessments of the effective-wave- 



length approximation, the validity of the Kirchoff 
model, and extensions of this work to optics having 
limiting apertures. One might also consider further 
simplifying physical and mathematical approximations 
in the diffraction theory used, as these can reduce 
computational resources needed for various applications 
[9]. 

Acknowledgments 

We thank A. C. Parr and D. G. Fischer for helpful 
discussions. Parallelized diffraction calculations em- 
ployed six Hewlett-Packard 9000 workstations coordi- 
nated with PVM (Parallel Virtual Machine) software. 
Version 3.3, developed at Oak Ridge National Labora- 
tory. 



6. References 

[1] W. R. Blevin, Diffraction Losses in Radiometry and Photometry, 
Metrologia 6, 39^4 (1970). 

[2] W. H. Steel, M. De, and J. A. Bell, Diffraction Corrections in 
Radiometry Opt. Soc. Am. 62, 1099-1103 (1972). 

[3] L. R Boivin, Diffraction corrections in radiometry: comparison 
of two different methods of calculation, Appl. Opt. 14, 2002- 
2009 (1975); Diffraction corrections in the radiometry of ex- 
tended sources, Appl. Opt. 15, 1204-1209 (1976); and Radio- 
metric errors caused by diffraction from circular apertures: edge 
effects, Appl. Opt. 16, 377-384 (1976). 

[4] L. R Boivin, Reduction of diffraction errors in radiometry by 
means of toothed apertures, Appl. Opt. 17, 3323-3328 (1978). 

[5] J. D. Purcell and M. J. Koomen, Coronagraph with Improved 
Scattered-Light Properties, Report of NRL Progress, 1962, U.S. 
GPO, Washington, DC (1962). 

[6] A. Rubinowicz, Die Beugungswelle in der Kirchoffschen Theo- 
rie der Beugungserscheinungen, Ann. Phys. 53, 257-278 (1917); 
K. Miyamoto and E. Wolf, Generalization of the Maggi- 
Rubinowicz Theory of the Boundary Diffraction Wave — Part I, 
J. Opt. Soc. Am. 52, 612-625 (1962), and Generalization of the 
Maggi-Rubinowicz Theory of the Boundary Diffraction Wave — 
Part II, J. Opt. Soc. Am. 52, 626-637 (1962); E. W Marchand 
and E. Wolf, Boundary Diffraction Wave in the Domain of die 
Rayleigh- Kirchoff Diffraction Theory, J. Opt. Soc. Am. 52, 761- 
767 (1962). 

[7] J. D. Jackson, Classical Electrodynamics, John Wiley and Sons, 
New York (1975) pp. 427. 

[8] T. Gravelsaeter and J. J. Stamnes, Diffraction by circular aper- 
tures. 1: Method of linear phase and amplitude approximation, 
Appl. Opt. 21, 364^3651 (1982). 



About the authors: Eric L Shirley and Raju U. Datla 
are members of the Optical Technology Division of the 
NIST Physics Laboratory The National Institute of 
Standards and Technology is an agency of the Technol- 
ogy Administration, U.S. Department of Commerce. 



753 



